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Abstract 

Calculations of L-shell x-ray absorption in transition metals are shown to be sensitive to screen- 
ing, both of the x-ray field and the photoelectron-core hole interaction. This screening is calculated 
using a generalization of the time dependent local density approximation and a projection onto 
a local atomic basis. The approach yields renormalized dipole-matrix elements which account for 
the observed deviations of the L^jLi branching ratio from the 2:1 value of independent electron 
theory. 
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Independent-electron theory is generally successful in describing near edge x-ray absorp- 
tion spectra (XAS) |TJ. However, it fails dramatically at the £2,3 edges in 3d transition metals 
0. |3|, f|, |^]. While the independent electron approximation predicts an L 3 /L 2 transition in- 
tensity "branching ratio" close to 2:1, the observed ratio (Fig. 1) varies considerably with 
respect to atomic number Z, and is closer to 1:1 for metals like Ti and V with nearly empty 
(i-bands || |7|, This puzzling behavior is thought to reflect many-body effects due to the 
Coulomb interaction, but despite many studies, its variation has never been quantitatively 
explained. We now show, however, that the observed branching-ratio depends crucially on 
dynamic screening of the x-ray field and the photoelectron-core hole interaction, and can 
be calculated using a simplified dynamic screening model and a generalization of the time- 
dependent local density approximation (TDLDA). Our approach makes use of a projection 
onto a local atomic basis and a real-space multiple-scattering (RSMS) formalism ||. This 
yields an efficient matrix formulation for extended systems, yet retains the simplicity of the 
TDLDA, and gives results in good agreement with experiment. 

The £2,3 XAS corresponds to transitions from the 2p\/2 and 2^3/2 levels to continuum 
s and d states. Several many body effects can be identified which contribute to a non- 
constant L3/L2 intensity branching ratio: i) Inelastic losses - these can be represented in 
terms of lifetime and self-energy effects in independent-electron calculations. The lifetimes 
are different for the L 2 and L 3 edges due to the Coster-Kronig mechanism JTD|, but this 
difference only increases the branching ratios, e.g., to about 3:1. ii) Dynamic core polarization 
- i.e., the creation of local fields which screen the external x-ray field. This polarization effect 



may be treated [y, |IT|] within the TDLDA by neglecting exchange terms, an approximation 
often referred to as the random phase approximation (RPA). This leads to a considerable 
reduction of the branching ratio, but does not account well for its variation with Z. iii) 
Screening of the photoelectron-core hole interaction - this effect, which we find to be crucial, 
can be addressed in terms of a frequency dependent exchange-correlation kernel f X c{^) in the 
TDLDA [Tj|, or by the analogous, non-local dynamically screened particle- hole interaction 
in the Bethe-Salpeter equation (BSE) |13|, [L4], [15], [I6| . The importance of dynamic screening 
of the core-hole is surprising, since it has been argued variously that corrections to the RPA 
are small |5|, [TT[] , or that an adiabatic kernel /° c is often adequate [[□], [T2| . 

The TDLDA [|5], [II], [T7| provides an efficient formalism for calculations of response func- 
tions, including corrections to the independent electron approximation, since it avoids the 
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FIG. 1: L3/L2 intensity branching ratio for the transition metal series from experiment (solid 
circles), and as calculated with different exchange correlation kernels: RPA (triangles), adiabatic 
f® c (diamonds); and the dynamic model of this work (see text) f xc (&) (squares). 

complications of non-locality in the time-dependent Hartree-Fock (TDHF) fl8|| , BSE, or 



configuration interaction approaches. The TDLDA was originally introduced for atoms |]TT 
but has since been extended to many other systems [|17] including band-structure formu- 
lations for transition metals ||. The TDLDA and TDHF equations are closely similar to 



the BSE [13|, |TJ], [Tj| [16], which provides a systematic, many-body framework based on the 
2-particle Green's function for treating core-hole screening. The crucial difference between 
the TDLDA and the BSE lies in the structure of the exchange-correlation kernel f xc (u); m 
addition the single-particle states are replaced by quasiparticle states which take the electron 
self-energy into account. Our approach (see below) is an extension of the TDLDA derived in 
part from the BSE. Like the BSE, our approach also uses a quasiparticle approximation for 
single-particle states, with inelastic losses approximated by a Hedin-Lundqvist electron-gas 
self energy ||19|| . Such losses are crucial in XAS and also point to the importance of dynamic 
screening. Our method simplifies the screening calculations by using a local basis. Local 
basis set methods have also been used to advantage in various related calculations |2C]. |2T| . 



Within the TDLDA |1 L] or TDHF |18] approximations, the photon cross-section (or XAS) 
<y(uo) can be expressed as an integral over the non-interacting response function xo(^^7 w ) 
and the screened x-ray field <j)(r, uj). For notational simplicity, it is convenient to regard the 
continuous coordinates r and r ' as vector or matrix indices, which may be suppressed unless 
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needed for clarity. Then the XAS can be expressed compactly as 



a(w) = -— ^(wJpmxoCwMw), (1) 



where 



Eq. (1) is equivalent to an analogous expression with replaced by the external x-ray field 
<fr ext and xo by the full response function \ Here fa are Fermi occupation numbers (1 or 
0), and the sums run over all one-electron eigenstates ipi(r) of the ground state Hamiltonian. 
The field <p(oo) consists of the external field <p ext = e ■ f (in the dipole approximation) plus 
an induced local field, which in matrix form is given by 

0M = e- 1 H0 e ^H, e(w) = l-tf(w)xV). (3) 

Here K(r,f,'uj) denotes the particle- hole interaction (or TDLDA kernel), which contains 
direct and exchange parts, i.e., 

K(f, Vu) = V(f,r') + f xc (f,f,'u), (4) 

and V = l/\r — r'\ is the Coulomb interaction. 

In this paper we consider several approximations for f xc (oo), which is generally a non- 
local, frequency dependent operator: i) The RPA (f xc = 0) - to the extent exchange can be 
neglected, the RPA is adequate ||. ii) Adiabatic TDLDA (f X c{0) = f xc ) - this static limit 
f xc (r, r ') = S(f—f')5v xc [p(f)]/Sp, is dependent on the local density and is obtained from the 
ground-state LDA exchange-correlation potential f xc [p]- hi) Dynamic TDLDA model - An 
LDA for the frequency dependence of f xc {u) was developed by Gross and Kohn [12]. At the 



large x-ray energies of interest here, this f xc (u) is strongly suppressed, and yields results close 
to the RPA ||. However, such results are clearly at odds with experiment for nearly empty 
(i-bands (Fig. 1). iv) Dynamic TDLDA/BSE model - Our aim here is to improve on i), ii) and 



iii) for L2,3 XAS, based partly on the BSE JT^, [14], [15], 0. In the BSE, the matrix elements 
(vc\f xc (uj)\v'c') depend on the dynamically screened Coulomb interaction W{uS) = e _1 (u;)l / , 
through an effective inverse dielectric matrix e _1 (co | ) ||15|| . However the actual dependence on 
uj ~ E c — E v is matrix element dependent, and depends on the effective dielectric response at 
the energy-transfer frequency, i.e., uj = uj + E c > — E v m E v > — E v . This behavior can be seen 
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explicitly in plasmon-pole models |13|]. For L 2 ,3 XAS, the most important occupied states 
v,v' are the 2pu 2 and 2p 3 / 2 levels, which are split by a moderate spin-orbit interaction A so , 
ranging from 5 eV for Sc to 20 eV for Cu. The matrix elements with zero energy transfer 
correspond to static screening; thus it is reasonable to set f xc {^) = f xc for v = v'. We also 
tried the unscreened, non-local TDHF exchange operator for v = v', which corresponds to 
an unscreened core-hole potential, but found it to be much too strong. For the off-diagonal 
elements (v ^ v'), however, we found that the unscreened exchange operator (i.e., the high 
frequency limit W = V) has only a small effect. This suggests that the effects of dynamic 
screening on off-diagonal terms at moderately high frequency ui = A so are also small and 
can be neglected. This behavior is in contrast to the case for optical absorption, where the 



adiabatic limit (cu = 0) is a good approximation and can be used for all matrix elements |L3 
Thus remarkably, we find that elaborate calculations of dynamical screening can be avoided 
for L 23 XAS by using a simplified dynamic model f xc {uj) — > f xc (£j), oj = E v — E v >, i.e., 
fxc{&) = fxc for = and f xc (&) =0, (uj — A so ) This defines our dynamic TDLDA/BSE 
model, which leads to reasonable agreement with experiment (Fig. 1). 

Next we briefly outline our calculations, which make use of the RSMS formalism (i.e., 
the real-space analog of the Koringa-Kohn Rostoker (KKR) band structure method) of our 
self-consistent, all-electron FEFF8 code [[5|. To begin we rewrite Eq. (1) as 

a{u) = M vL {u)p L , L ,{E)M vL ,{u), (5) 

C v,LL> 

where E = u + E v — Ep is the photoelectron energy. The screening of both the x-ray 
field and the photoelectron-core hole interaction are included implicitly in the renormalized 
dipole matrix elements |fL8| , M vL {uj) = (R L \<p\v), where L = (k, m) denotes a relativistic 



angular momentum basis. The quantities pL,v{E) are matrix elements of the unoccupied 
one-electron density matrix, 

p(f, r ', E) = Mr)€(f')S(E - E c ), 

c 

= Y R L(r)RL'(r')pL,u(E), 

L,L' 

p L> v(E) = 5 L:L , + xl,l>(E). (6) 

Here Rl{^ , E) are normalized scattering states calculated with the absorbing atom potential, 
and xl,l'{E) contains the fine structure in the XAS due to scattering by the environment 
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||. Note that by replacing <fi with <p ext in Eq. (5), the screened dipole matrix elements M v l 
become bare dipole matrix elements M v l = (Rl\£ • and one recovers the independent 
electron formula, equivalent to Fermi's Golden Rule. Since the strength of the XAS is a 
measure of the screening response, the independent-electron approximation should become 
increasingly valid away from the edge region. 

The second key approximation in our approach is the use of a local basis for calculations 
of Xo an d M v l. This is done starting from an expression in terms of a Kramers-Kronig (KK) 
transform over the density matrix, 

dE 



Y,r v (r)Mr') / —p{?,r!E) (7) 



1 



+ 



lu- E + E v + i5 u + E - E v + i5 



Once xo is known, Eq. (|3]) could be solved iteratively in real space to obtain 4>(r) (TT|. How- 
ever, this procedure is computationally expensive for extended systems, since it involves KK 
transforms for many (r, r' uj) points. To simplify this calculations, we make the reasonable 
assumption |TT| that the induced charge p md = xo(^)0 that screens the x-ray field, is local 
and arises largely from a few significant orbitals on the absorbing atom. This is convenient, 
since our formulation only needs the screened field (j>(r,u) at short distances to calculate 
the deep-core transition matrix M v l{oj). Thus to approximate 0, we introduce the atomic 
projection operator P = J2 n \i ) n)('4 ) n\, which projects a given function onto a local basis set 
of atomic-like orbitals on the central atom. Then the density matrix can be approximated 
by its local contribution p w p loc = PpP. These approximations can be systematically 
improved by including a more complete set. Thus 

dritt'u) = £C(0C(Ox™VH (8) 

vnn' 

ioc i , .\ k r ^{n\R L )p L:V (R L ,\n') 



— Y dE- 

7T j~t, Je f uj — E + E v + id 

where k = \Jl{uj + E v ). Note that the localized part of xo does not require a KK transform 
at each point, since the localized part of the photoelectron wave function can be separated 
into energy and position dependent parts. Moreover, the overlap matrices {v\Rl) decay 
rapidly with energy, so the KK transform converges well. This approximation then leads 
to a fast matrix formulation for M v l- From Eq. (|3|), we obtain (summation over repeated 



indices being implicit) 

(j){r,u) « (f) ext (r, lu) + ^2 J dr' K(r,r,' uj) 

v'n'n" 

x </(r')^'(r') xl'yn''^)^'"" ) (9) 

where M vn = (^ n |^»|i;) is calculated by integrating Eq. @ over the core and basis set 
functions, 

M m (w) = M vn + K m y n > Xv'n'^'n"^) AW'( W )> 

M m = (n\(j) ext \v), K m>v ' n i = (vn\K\v'n'). (10) 

These equations can readily be solved by matrix inversion. Finally on integrating Eq. (|9p 
over the core- and final state-wave functions, we get 

M vL {i0) = M vL (iu) +K vL>v , n ,xin>,v'n»M v > n >>, 

K vW = (vR L \K{u)\v'n'), (11) 

where Rl denotes the scattering-state Rl{^ — E v ). The matrix form in Eq. ( pf) is very 
efficient, and has been implemented using an extension of our FEFF8 code. This extension 
is straightforward, since only the dipole matrix elements need be modified to incorporate 
screening. Due to the local form of f xc (uj), the contributions to K v iy n '(u) satisfy the same 
selection rules and can be calculated using standard formulas for Coulomb interaction matrix 
elements p2| . 

Typical results near the beginning and end of the transition metal series are presented 
for Ti and Ni in Fig. 2. The dramatic differences reflect differences in the response between 
nearly empty and nearly filled <i-bands, and are strongly dependent on the form of fxci^)- 
For all calculations we used theoretical atomic core- hole life-times |TI|. We did not add 



additional broadening to correct for experimental resolution, though this would give slightly 
better agreement with experiment (Fig. 1). Our results for the RPA agree well with those 
of Ref . H , which validates our local screening approximation. Note that the RPA is only a 
good approximation for nearly empty <i-bands, while the adiabatic /° c is appropriate only 
for nearly filled ones. However, our dynamic model f xc (uj) is clearly satisfactory for the 
entire series. 

Since screening redistributes the oscillator strength between the L 2 and L 3 edges, the 
importance of these many-body corrections appears to cast doubt on the accuracy of results 
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FIG. 2: L3.2 edge XAS vs energy with respect to the Fermi level, for a) Ti (upper figure) and b) 
Ni (lower), as calculated with different screening models: RPA (solid), static fxc° (dots), and the 
dynamic model of this work f xc (ui) (dashes). 

obtained from the XAS "sum- rules" |J . These sum-rules allow one to determine various spin- 
and orbital-moments from linear combinations of the L2 and L3 XAS [[23] . However, one can 
now correct these procedures for local-field effects with our approach, e.g., by substituting 
the screened XAS cross-section, in place of the one-electron result in the analysis. 

In summary, we have found that dynamic screening of the photoelectron-core hole inter- 
action, which gives rise to a frequency dependent exchange-correlation operator f xc (w), is 
crucial in calculations of transition metal £2,3 spectra. However, we have found a dynamic 
model based on the TDLDA and BSE, which accounts well for the frequency and matrix 
element dependence, by neglecting off diagonal, high frequency screening terms. With this 
model, we have developed an efficient approach for including screening in deep-core XAS, 



S 



based on calculations of screened dipole matrix elements. Our approach goes beyond the 
conventional TDLDA, and is similar in some respects to a screened TDHF approximation 
16| . Moreover, the approach yields good agreement with experiment for the L 2i 3 XAS of 3d 



transition metals, without the complexity of full dynamic-screening calculations. 
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